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We describe a proposal for fast electron spin initialization in a negatively charged quantum dot 
coupled to a microcavity without the need for a strong magnetic field. We employ two-photon 
excitation to access trion states that are spin forbidden by one-photon excitation. Our simulation 
shows a maximum initialization speed of 1.3 GHz and maximum fidelity of 99.7% with realistic 
system parameters. 



Recent demonstrations of cavity quantum electrody- 
namics (QED) effects with semiconductor quantum dots 
(QDs) coupled to microcavities [1, 2, 3, 4, 5] show that 
these systems are robust and scalable platforms for quan- 
tum information science. The QD-cavity QED exper- 
iments performed so far treat each QD as a two-level 
quantum system consisting of a ground state and the 
single exciton excited state. However, multilevel quan- 
tum systems are necessary for fast, complex, and high 
fidelity quantum information processing [6]. Each of 
these systems should contain two stable ground states 
between which the transition is dipole forbidden and one 
or more excited states that serve as a passage for pop- 
ulation transfer between the ground states (A system) 
[7, 8]. Initialization, coherent manipulation, and readout 
of these A systems are essential for quantum informa- 
tion processing including controlled phase gate [9, 10], 
quantum repeaters and networks [11], and remote entan- 
glement distribution or creation [12]. Therefore, realizing 
a multilevel system in a QD coupled to a cavity is crucial 
for semiconductor cavity QED implementation of quan- 
tum computers. 

One way to realize a A system in a negatively charged 
QD is by employing the Zeeman-split electron spin states 
as the ground states and a trion state as the excited state 
[6]. Recent experiments with QDs (not coupled to a cav- 
ity) have demonstrated spin initialization with magnetic 
field along the QD growth direction (Faraday geometry) 
[13] as well as with field perpendicular to the QD growth 
direction (Voigt geometry) [7, 14, 15]. Initialization in 
the Faraday geometry requires weak magnetic field and 
can achieve very high fidelity, but the initialization is slow 
because the process relies on random spin-flip Raman 
scattering. On the other hand, initialization in the Voigt 
geometry is fast, but a strong magnetic field is required 
to mix the spin states and split the resulting cigenstates 
for high fidelity. The high cost and large space required 
to achieve a strong magnetic field present a significant 
drawback of this initialization method. 

In this Letter, we propose a fast and high fidelity spin 
initialization method for a negatively charged QD cou- 
pled to a microcavity. Our method also works for a pos- 
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FIG. 1: (color online) (a) Quantized energy levels of a neutral 
QD. The potential in the xy plane (perpendicular to the QD 
growth axis) is modeled as a two-dimensional harmonic os- 
cillator. Upon two-photon absorption, an electron is excited 
from the valence band s orbital (\g)) to the conduction band p 
orbital (|e)). (b) Energy level diagram of relevant electron and 
trion states in a negatively charged QD. Electrons (holes) are 
shown in solid (opened) circles. Subscripts indicate orbital 
labels. The ground levels are | f s ) = |1) and ] J. 8 ) = 2). 
The first excited state is the trion state | tslsfTs) = |5}- The 
second excited states are trion states excited by two-photon 
processes | TsTp1%) = 1 4) and | J. s tplTs) = |3). The straight 
arrows indicate two-photon excitation, the curly arrows indi- 
cate various decays of excited states, and the curved arrow 
indicates spin flip. 



itively charged QD [16], but we will focus on the nega- 
tively charged one as an example in this paper. The ma- 
jor advantage of our method is the absence of a strong 
magnetic field, because there is no need to mix the elec- 
tron or the hole spin states. Instead, we employ two- 
photon process to excite transitions that are otherwise 
spin forbidden by one-photon excitation. A doubly reso- 
nant cavity is used to enhance both the excitation and the 
spontaneous emission rates. We show that our method 
can achieve a maximum initialization speed of 1.3 GHz 
and fidelity of 99.7% with realistic system parameters. 

We consider a strongly confined QD and neglect the 
mixing between the light hole and heavy hole bands be- 
cause of their large energy splitting. We model the QD 
potential as a finite well along its growth axis (z axis) and 
a two-dimensional oscillator perpendicular to its growth 
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axis (xy plane) [17, 18]. The effective mass for electron 
(hole) is to* {m* h ) and the oscillator frequency for the 
conduction (valence) band is uj e (u>h)- The harmonic os- 
cillator quantized levels are labeled as s orbital, p orbital, 
etc. (Fig. la). Because the QD is much more confined 
along the z axis than in the xy plane, the first several ex- 
cited states have wavefunctions occupying different states 
in the xy plane but the same ground state along z axis 
[19]. 

Fig. lb shows the relevant QD energy levels for our 
initialization method. The ground levels arc the electron 
spin up | |s) and spin down | J, s ) in the s orbital. The 
application of a weak magnetic field of 0.2 T reduces 
hyperfine induced spin-flip rate 712 but does not mix the 
two spin states [20, 21]. The first excited state is the 
trion state | tslstls) with the electrons in the conduction 
band s orbital and the hole in the valence band s orbital. 
The second excited states are trion states generated by 
two-photon excitation | tstpfls) and | Istpfl's) in which 
one of the electrons is in the conduction band p orbital. 

The spin initialization to state | "f s ) is performed as 
follows: Two-photon excitation with a + circularly polar- 
ized light transfers population from | "f s ) to | Tstpfl's) and 
from I J, s ) to I Istpfl's) with a rate f2. The electron in the 
conduction band p orbital then relaxes to the conduction 
band s orbital by a phonon assisted process with a typical 
timescale of 10-30 ps [19]. Due to Pauli exclusion prin- 
ciple, the resulting state is | Ts-UfTs)- The population in 
this state will preferentially decay to the state | ts) with 
rate T = 715 because the decay to state | | s ) is spin for- 
bidden. Eventually, the population accumulates in state 
t s ). Two-photon excitation cannot directly create an 
electron and a hole in the conduction and valence band 
s-orbitals as these two orbitals have different parity. For 
this reason, we use p-orbital for intermediate states. 

In order to achieve fast initialization, a doubly reso- 
nant cavity is needed to enhance the QD spontaneous 
decay rate T and the two-photon excitation rate il. We 
consider InAs/GaAs QDs with s-s transition wavelength 
of 930 nm. Since the p orbital is approximately 50 meV 
above the s orbital in the conduction band [19], this im- 
plies that the s-p transition wavelength is approximately 
900 nm. Therefore, we need a doubly resonant cavity 
at 930 nm and 1800 nm to enhance the spontaneous de- 
cay rate and the two-photon excitation rate respectively. 
Both resonant modes should be polarization degenerate, 
because all the selection rules described in this Letter 
require circularly polarized light. 

We follow Ref. [22] to estimate the two-photon exci- 
tation rate of a QD inside a microcavity. When both 
photons are of the same frequency, the cavity-enhanced 
two-photon absorption rate is T^ PA = 2ir\Q\ 2 5(u>d — 2cj), 
where the effective Rabi rate f2 is 

~ ^ A gk 



where 77, P, Q, n, and V are the light coupling efficiency 
into the cavity, excitation power (measured outside the 
cavity), cavity quality factor, refractive index of material 
that the cavity is composed of, and cavity mode volume, 
respectively. tu c , uj^, and to are the cavity resonance fre- 
quency, the transition frequency between two QD levels 
coupled by the two-photon process, and the laser fre- 
quency, respectively. d gk (d ke ) is the dipole moment 
between the ground (excited) state and the intermedi- 
ate state k. A gk — u> k — u> g — u> is the laser detuning 
from the transition between the QD ground state and 
the intermediate state. The delta function 5{u>d — 2w) 
represents energy conservation. The factor <f> describes 
the spectral mismatch between the laser and the cavity, 
cj>(ui) = 1+4Q2 "f,/^ e _ 1) 2 ■ The factor ip gk (ip ke ) character- 
izes the reduction in Rabi rate due to position mismatch 
between the QD and the electric field maximum within 
the cavity, ^ gk = ^ggL (^f), where E(r), E(r M ), 
and e are the electric field at r, electric field maximum, 
and electric field polarization at the QD location, respec- 
tively. 

For the transition from valence band s orbital to con- 
duction band p orbital, the dominant intermediate states 
are the conduction band s orbital \k\) — \s c ) and the 
valence band p orbital |fc 2 ) = \p v ) (Fig. la). In our 
analysis, we assume |d gfcl | = |d efc2 | = e\v cv \ {\v cv \ = 0.6 
nm [23]), which is the transition moment between the 
valence band and the conduction band integrated over a 
unit cell. We also calculate \d ekl \ — el e and |d 9 fe 2 | = elh, 
where e is the electronic charge and l e (lh) is the oscillat- 
ing length of the electron (hole) with l e — yjhj {2m*Lo e ) 
(l h = y/h/(2m* h uj h )). With rj = 2%, Q = 5000, 
V = (A/n) 3 , <f> = 1, ip = 1, and P = 50 /iW, we 
calculate H,/2tt — 4.5 GHz. We note that we use very 
conservative number of 2% for in-coupling efficiency into 
the cavity, as in our previous experiments on controlling 
cavity reflectivity with a single QD [1], but significantly 
higher coupling efficiency exceeding 80% is possible with 
fiber tapers or cavity- waveguide couplers [24, 25]. There- 
fore, much higher £1 is possible with a smaller excitation 
power. 

The fidelity and speed of our spin initialization method 
are analyzed using method based on the master equation 
[26]. Here we analyze the case when er + light is applied. 
The time evolution of the density matrix p of the QD is 
given by 
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with the Hamiltonian 

n = m{\ TsUsKT* I + 1 UplMU. I) + h - c - 

The Zeeman splitting between the states | f s ) and | J. s ) is 
approximately 5.4 /xeV for an electronic g-factor of 0.46 
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[15] and a magnetic field of 0.2 T. This splitting is much 
smaller than Ml and hT, therefore we neglect it in our 



analysis. Cj(p) = DjpD^ 
.D's are the collapse operators: 
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hpD'jDj, where 
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D 2 = Vl25\ ls)(hlsts |, 

D 3 = \/735| TsisfTs) (IsTpfTs |, and 

D4 = v/TiFI T.i.ft.XUpfta I- 

Numerical values of the rates used for the simulation are 
725 /2tt = 100 kHz [13], 7i 2 /2tt = 10 kHz, and 735/2^ = 
745/2-71 = 8 GHz corresponding to a relaxation time of 
20 ps (Fig. lb). The QD spontaneous emission rate T 
is Purcell enhanced and its numerical value ranges from 
5 to 20 GHz, corresponding to a Purcell factor ranging 
from 50 to 200. 

The fidelity and initialization speed have strong de- 
pendence on f2 and T. Experimentally, f2 is tuned by 
changing the excitation light intensity and by employing 
cavity resonance at two-photon excitation frequency; T is 
tuned by placing the QD in nanocavities of various qual- 
ity factors. To find the fidelity, the steady-state density 
matrix is solved. The fidelity is defined as the popula- 
tion in state | | s ), namely p\\. To find the initialization 
speed, the density operator is evolved in time from the 
initial condition of pu(t = 0) = P22(t = 0) = 1/2, where 
P22 is the population in the state | | s ). The initialization 
time is defined as the time required for pn(t) to reach 
the value of 1 — 1/e. 

Fig. 2a shows the initialization fidelity as a function of 
the two-photon effective Rabi rate f2 and the one-photon 
spontaneous emission rate T. For fixed T, increasing £1 
depletes population from state | | s ) more efficiently and 
therefore leads to increasing population in state | | s ) and 
thus higher fidelity. However, while depleting population 
from state | J. s ), two-photon excitation also depletes pop- 
ulation from state | j s ) (Fig. lb). Therefore, there are 
two competing processes: one depletes the population 
in I | s ) with a rate of f2 and the other populates | 1 S ) 
with a rate of T. Hence, increasing f2 indefinitely, for a 
fixed T depletes | | s ) and ultimately decreases the fidelity 
of initialization (Fig. 2a, b). Such reduction in fidelity is 
more prominent for lower T. On the other hand, for fixed 
fi, increasing T allows faster population transfer from 
state I TsJ-sHs) to state | f s ), leading to higher fidelity. 
The maximum fidelity obtained is 99.7%. This fidelity is 
comparable to the fidelity of initialization achieved in the 
Faraday geometry (> 99.8%) [13] at a magnetic field of 
0.2 T and in the Voigt geometry (ss 99.7%) at a magnetic 
field of 1 T [15]. It should be noted that the spin flip rate 
712 also limits the fidelity. For a spin flip rate increased 
from 7i 2 /2tt = 10 kHz (as in Fig. 2) to 100 kHz and 1 
MHz, the maximum fidelity decreases to 99%, and 98% 
respectively. 

Fig. 3a shows the time needed to perform spin initial- 
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FIG. 2: (color online) (a) Fidelity of electron spin initializa- 
tion as a function of two-photon excitation rate Q/2ir and 
spontaneous emission rate T/2ir. (b) Slices through plot (a) 
for specific T/2n values, indicating maximum attainable fi- 
delity of 99.7%. 
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FIG. 3: (color online) (a) Time needed to perform electron 
spin initialization as a function of and F. (b) Initialization 
time as a function of the ratio of SI and F for different values 
of r. 
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ization as a function of f2 and T. As shown in Fig. 2b, 
an increase in fl/T results in decrease in the fidelity . In 
fact, for il/T > 0.2, the fidelity never reaches the value 
of 1 — 1/e and therefore initialization time is arbitrarily 
long. For this reason, we study only the range fi/T < 0.2 
in Fig. 3b. It should be emphasized that the value 0.2 is 
obtained for the parameters used in the simulation, and 
a different value will be obtained for a different set of 
parameters. As shown in Fig. 3b, when f2/T < 0.2, the 
initialization speed increases with increasing f2/r. Al- 
though increasing T in general increases the speed, this 
speed increase is not significant when T 735 = 745 be- 
cause the phonon assisted relaxation rates from | TsTpfN) 
and I IsTplTs) to I TsisfN) become the limiting factor. 
The minimum initialization time calculated is 120 ps, 
giving the maximum speed of 1.3 GHz. This speed of 
initialization is much higher than the maximum speed 
obtained in the Faraday geometry (100 kHz) [13] and 
the Voigt geometry (144 MHz) [15]. It should be noted 
that the speed in Voigt geometry can potentially be en- 
hanced by using a cavity, but that is not true for Faraday 
geometry. Unlike fidelity, a change in spin-flip rate 712 
does not affect the speed of initialization, as 712 << ^, T. 

It is important to note the trade off between high speed 
and high fidelity. Fig. 3 shows that high two-photon 
excitation rate f2 is needed to achieve high speed, but 
subsequently the fidelity decreases at high £1 (Fig. 2). 
Therefore, the speed of initialization corresponding to the 
maximum fidelity of 99.7% is lower (« 50 MHz). Sim- 
ilarly, at the highest initialization rate, the fidelity of 
only 80% can be achieved. However, we can still achieve 
moderate fidelity with moderate speed: for example with 
H/2tt = 0.5 GHz and T/2ir = 10 GHz we can achieve a 
fidelity of 97.3% and a speed of 320 MHz. 

In conclusion, we have described a new method for 
initializing electron spin in a negatively charged QD us- 
ing cavity enhanced two-photon excitation. The highest 
speed of 1.3 GHz and the maximum fidelity of 99.7% can 
be achieved with weak magnetic field of approximately 
0.2 T and realistic system parameters. Although it is not 
possible to achieve both the maximum speed and fidelity 
simultaneously, with H,/2tt — 0.5 GHz and r/27r = 10 
GHz, one can achieve a fidelity of 97.3% together with 
a speed of 320 MHz. The use of weak magnetic field re- 
duces the cost and the space requirements for a quantum 
network, as magnetic field of 0.2 T can easily be achieved 
with a small electromagnet, as opposed to using a large 
superconducting magnet needed for strong magnetic field 
in other spin initialization techniques. 
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